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Abstract 

We evaluate the low-temperature cutoff for quantum interference 1/t<^ induced in a d-wave 
superconductor by the diffusion enhanced quasiparticle interactions in the presence of disorder. We 
carry out our analysis in the framework of the non-linear cr-model which allows a direct calculation 
of 1/t<^, as the mass of the transverse modes of the theory. Only the triplet amplitude in the 
particle- hole channel and the Cooper amplitude with is pairing symmetry contribute to We 
discuss the possible relevance of our results to the present disagreement between thermal transport 
data in cuprates and the localization theory for d-wave quasiparticles. 

PACS numbers: 74.25.Fy, 72.15.Rn 
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I. INTRODUCTION 



During the last ten-fifteen years overwhelming experimental evidences for a li— wave pair- 
ing symmetry in high-T^ superconducting cuprates have been collected^. Such an uncon- 
ventional symmetry plays a crucial role in determining the physics of the superconduct- 
ing phase. In particular the presence of four nodes in the (i— wave superconducting gap 
Ak = A(cosA;^. — cosky)/2 and, consequently, the existence of gapless quasiparticles (QP) 
excitations down to zero energy strongly affects the low temperature transport of these sys- 
tems. The energy spectrum of the QP in the proximity of the d— wave gap four nodes has 
a Dirac cone shape 



= + ^ ^v^l + vlkl (1) 

where is the energy dispersion of the non interacting electrons, while vp and f a represent 
respectively the Fermi velocity and the slope of the superconducting gap Ak at the nodes. 
A weak amount of disorder originates^ a finite quasiparticle elastic lifetime tq and a finite 
density of states at the Fermi level A^o = l/(7r^f_F'i^A''o) ln(po''o)) Po being a cutoff of order 
A . In the low temperature and frequency regime the conductivities of the system (both 
the electrical a and the thermal k) do not depend on the amount of disorder but only on 
the bare QP spectrum parameters vp and f a- This leads to the so-called universality of the 
low energy values of the conductivities in c?— wave superconductors. Indeed, by allowing the 
impurity-scattering between all the four nodes of the d— wave spectrum and neglecting the 
vertex corrections, the following "universal" expressions are found^ 

hm ain^T) = 4^, hm = ^>1±A. (2) 

The inclusion of vertex corrections modifies cr by a factor depending on disorder and Fermi- 
Liquid parameters, while it leaves k totally unchanged^. 

However, in cuprate superconductors this generic result could be strongly affected by 
localization corrections because these materials are quasi-bidimensional and quantum inter- 
ference effects are usually relevant in low-dimensional systems^. 

/^From a theoretical point of view the evaluation of the localization corrections to the 
conductivities in a (i— wave superconductor is a really complex topic. Indeed, when consid- 
ering a (i— wave BCS Hamiltonian in the presence of disorder and neglecting interactions 



among quasiparticles, the theory for the quantum interference predicts a multitude of dif- 
ferent regimes and crossovers depending on the specific symmetries of the underlying Fermi 
surface and of the disorde r^i^i^i'^'i^i^i-^" . Anyway, in the generic case of non-magnetic disorder 
connecting all nodal points the theoretical calculations indicate the occurrence of an Ander- 
son localization of the c/— wave QP excitations^. The theory also predicts precursor effects, 
the so-called weak localization corrections, which are logarithmic in temperature. 

The analysis of experimental findings for the thermal transport in the cuprates does not 
appear to fit the above theoretical scenario. Indeed when comparing these predictions with 
the experimental measurements a convincing agreement is found neither with the "universal 
behavior" nor with the complete localization. It is worth stressing that thermal conductivity 
is conveniently used in order to test the theoretical results, because, differently from charge 
conductivity, it is not affected by the vertex corrections allowing for a direct comparison 
between the expression of n{T) in Eq. Q and the experimental data. 

Thermal conductivity measurements have been performed on various cuprates, down to 
temperatures of the order of 100 mK. In particular the residual thermal electronic conduc- 
tivity Kres/T (= limT_>oK/T), extrapolated from measurements carried out on optimally 
doped YBa2CusOe+^ {Y B CO (123))^^'^^ and Eis^rsCaCnaOg (BSCCO)^^- seems to be 
in agreement with the universal result (Eq. (0) with vp and va extracted from ARPES 
measurementsi^ii^) . The agreement is particularly good for BSCCO, leaving no room for 
sizeable weak localization corrections. The situation is partially different for YBCO{123), 
where the inclusion of sizeable weak localization corrections would provide not unplausible 
larger ratios vf/va than the estimate of Ref. However (or moreover) thermal mag- 
netoconductance measurements^ do not give any indications of the presence of detectable 
weak localization corrections which should appear as a positive crossing to a InH con- 
tribution at very small fields. This result is quite generic for all materials at different 
dopings. On the other hand, measurements performed on YBa2Cu40s {YBC0{12A)) un- 
derdoped and Pr2-xCexCuOi (PCCO) optimally doped provide a residual conductivity 
Kres/T much Smaller than the universal value. Such a result, compatible with Hres/T = 
could be ascribed to the localization of the quasi-particles responsible of the low tempera- 
ture energy transport. Recently measurements of Hres/T at various doping became available 
in La2-xSrxCu04^ (LSCO, Ref. Il6[ ) and in YBCO{123y^ and made the theoretical sce- 
nario even more involved. In particular a significant doping dependence of Kres/T has been 
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observed in these cuprates. Such a dependence is hardly compatible with the universal 
behavior, particularly in the LSCO compounds, since it would imply a strong variation on 
doping of vp and f a in partial contrast with ARPES measurementsi^ and the tempera- 
ture dependence of the penetration depth. Moreover ^res/T turns out to be finite and no 
localization is found for all dopings in the superconducting phase. 

Concerning the absence of weak localization contributions a main issue is the estimate of 
the temperature below which quantum interference effects start to be visible before local- 
ization occurs. In the metallic phase this issue involves the determination of the dephasing 
time T^p, i.e. the time scale above which coherence is destroyed by inelastic processes. The 
energy scale represents a cutoff for quantum interference and plays a crucial role in 
determining quantitatively the magnitude of the localization effects at low temperature^. 

In principle, the situation becomes more involved when considering a (i— wave supercon- 
ductor. Due to the extra symmetry of the Cooper pairing, quantum interference determines 
a relevant correction to the density of states which is logarithmic in energy since it derive 
from integration on diffusive modes with energy as a cutofi^*^^. This fact reflects in the 
conductivity with the onset of an other source of logarithmic corrections coming from the 
thermal average of the density of states contribution. These corrections no longer depend 
logarithmically on the dephasing time (oc — Z?2(r^/ro)), but directly on T similarly to the 
interaction contributions^ (oc — Zn(l/Tro)). As a consequence one can naturally individuate 
two contributions to the conductivity corrections, coming out respectively from the modes 
damped by t^^ (as in the normal metal) and from those damped by the energy 

that can be viewed, roughly speaking, as the renormalization of the diffusion coefficient 
and of the density of states {5as ~ 6DNq + D6Nq). Indeed half of the contribution to the 
Renormalization Group (RG) equation^ for ag comes from the first term and half from the 
second term. 

However it can be easily shown that in the range of temperature from lOOmK to IK, to 
which we are interested in, the second kind of corrections can be neglected!-^, because of the 
rather small value of I/Ttq which is usually much less than t^/tq. We can therefore estimate 
the localization correction to the spin conductivity as of a (i— wave superconductor, (and 
consequently to the thermal conductivity which is related to by a generalized Wiedmann- 
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Franz law^), considering only the expression for (— given by 



^O's I _ Stores _ ^ ly^^f 
as /^-i l^res 2 To 

where t = Xjli^^Os and Og = {vp + vj^)/{Tr'^vpVj\). 

Note that the right term of the Eq. 01 differs for a factor 1/2 from the analogous RG 
expression^ since, according to the above discussion, only in half of RG contribution the 
relevant time region for the interference ranges from Tq to r<^, while in the other half it 
ranges from tq to 1/e ~ T. 

In normal metals the most effective mechanism of dephasing is provided by the interac- 
tions among electrons^S. Indeed interactions in presence of disorder have a twofold effect: on 
one side they generate corrections which can compete with those due to pure quantum inter- 
ference, on the other they provide an intrinsic dephasing time T^p which limits the quantum 
interference processes for T 7^ 0. 

A systematic evaluation of interaction corrections to the conductivities in the d-wave 



y (see alsoQ. 



Since these corrections 



superconducting phase was carried out in Refs. 
are proportional to the disorder induced density of states A'o, they should be negligible for 
clean enough systems and for temperatures of the order of 100 mK, at which the thermal 
transport measurements are presently available^. 

The aim of this work is to extend the analysis of Ref . and to derive an explicit expression 
for the dephasing time induced in the superconducting phase by the interactions among 
the rf— wave quasiparticles. In particular we will determine 1/t^ as the mass of the diffuson 
and the Cooperon (i.e., the two particle propagators in the particle-hole and particle-particle 
channels) in analogy with Refs. 2ll)22 22 where 1/r^ was evaluated in the metallic phase. 
Here we will perform a one-loop analysis of the interaction terms, within the non-linear a- 
model which is described in detail elsewhere^i^. We find that the most relevant contribution 
to 1/r^ comes from the triplet channel. Indeed the singlet channel contribution is ruled out 
by the lack of particle conservation in the superconducting stat o^i^^ . On the other hand the 
residual interaction in the Cooper channel, which would lead to a d + is instability when 
it is attractive, provides a subleading contribution to 1/r^ when it is repulsive, as we shall 
assume here. As a final outcome the following estimate for ratio t^/tq, entering in Eq. El 
will be derived 

^^^(Aro). (4) 

To J 



where ep is the Fermi energy and the proportionahty factor is of order one. In agreement 
with the resuhs for the metalhc phase is proportional to the inverse of temperature 1/T. 
Moreover the hnear dependence of the ratio Ti^/to on the elastic scattering time suggests that 
quantum interference corrections at finite temperature would be larger in cleaner systems, 
differently from what one could expect intuitively^^. 

By inserting in Eq. Q the experimental determinations of vp, va, A and ep, we will 
estimate the value for different cuprates, in order to predict from Eq.Q the temperature 
below which localization corrections should become relevant. In particular in YBCO{123) 
and BSCCO, optimally doped, we will find that the corrections to the universal results 
should become relevant, of order of 20 -ir 30 per cent, for temperatures of the order of 0.1 -^ 1 
K, in disagreement with the experimental results. Actually, it is not clear why localization of 
the quasiparticles does not occur in YBCO{123), while an explanation can be found for the 
bismuth cuprate BSCCO, as we will explicitely point out in Sec. IV. Indeed a rough estimate 
of the transport times, from electrical conductivity at higher temperature, suggests the 
existence in this last compound of other mechanisms for decoherence and dephasing, more 
efficient than the mechanism of interaction in a disordered environment we have considered 
here. 

Finally, in sec. V we will present the concluding remarks on the problems posed by our 
findings. 

II. FIELD THEORY FOR DISORDER IN THE SUPERCONDUCTING PHASE 

In this section we summarize the field theory approach and introduce the non-linear a- 
model (NLcrM), that will be used in the following to evaluate the dephasing time I/t^^ in the 
(i— wave superconducting phase. It is worth noting that in the superconducting phase the 
evaluation of the dephasing time using the NLcrM is simpler and more direct than using the 
standard perturbation theory, because of the presence of both the normal and the anomalous 
Green functions with specific k dependences. 

For the sake of clarity, we start discussing in this section the NLaM in the absence of 
quasiparticle interactions, which will be included in a second step. In order to describe the 
properties of the disordered system within a path integral formulation, we introduce the 
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spinorial representatio: 
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= ^ I '^'^ ^ I , ^{r) = [C^^. (5) 
V ^ \ iayc{r) 

\&(r) is an Nambu spinor, while c(r) and c(r) are Grassmann variables with components 
c{r)(j^p^a and c{r)(j^p^a, where a refers to the spin, and p = ± is the index of positive 
and negative {—fi) frequency components. Later p will be extended to label the positive 
and negative Matsubara frequencies Um = TcTm, where m is an odd integer, a = l...n is the 
replica index introduced by the replica trick used to average the action over disorder, and 
C is the charge conjugation matrix 

C = iayTi. (6) 

Here and thereafter, the Pauli matrices ai{i = x,y,z) will act on the spin components, 
Si{i = 1,2,3) on the frequency components, and Ti{i = 1,2,3) on the components of the 
Nambu spinors Cj and q. 

In this representation the BCS action, describing the system in the superconducting phase 
in the absence of disorder, reads 

5o = J] (a + i^kT2Si - ins3) (7) 

k 

where the term —iQs^ has been introduced in order to define retarded and advanced Green 
functions, and the presence of the Pauli matrix si in the second term is due to the fact that 
the Cooper pairing couples states with opposite energies. 

As in the standard theory of localization in the normal phase^^, disorder is introduced 
via an impurity scattering potential V{r) with local Gaussian distribution and zero mean 
value (the overline indicates the impurity spatial average) 



V{r) = , V{r)V{r') = u'6{r - r ) 

In the following we will write u"^ = 1/ {2nNQTo) as set by the saddle point solution (see 
below). 

The average over disorder distribution is performed by using the replica trick method. A 
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four field term is generated which couples fields with different replica indices 



^^^. = -^/^^(^(l/)^(l/))' (8) 
This contribution is then decoupled with a standard Hubbard-Stratonovich transforma- 
tion. Consequently we introduce an hermitian bosonic field Q, Q being a matrix in the 
(cr, r, p, a)-space. The decoupled action has the following expression 

Simp = ^ {/ I^NoTtQ' - Q^)] I 

where Q"^) represents a scalar product. 

The NLcrM can now be derived, after an integration over the Grassman variables, by 
expanding the effective action around the saddle point Qsp = l/(2ro)s3 and taking into 
account only the low energy transverse fluctuations. As a consequence, the fluctuating fields 
Q satisfy the conditions = 1, TiQ = and can be written in the form 

Q = U-'QspU, (9) 

where f/ is a unitary matrix with suitable symmetries. We also rescale Q 1/2tqQ in order 
to express the action in term of an adimensional field. 
The outcome of this procedure is the following NLcrM 



SsM = ^ 



!^as j dxTr{VQY (10) 
8Non J dxTris^Q) 



where 

nV ^ (E2 + (2ro)-2)2 ' vpv^ ^ ^ 

is the quasi-particle conductance in the Born approximation which can be interpreted as the 
spin conductance of the systen>^. The appearance of ag (instead of the charge conductivity, 
as in the metallic case) in front of the term Tr(V(5)^ in Eq. (|l(Jj) individuates the inverse 
spin conductivity as the natural expansion parameter of the theory for the superconducting 
phase. 

The expressions for the diffuson and the Cooperon (which are the massless modes of the 
theory showing diffusive poles) can be calculated by expanding the NLcrM action Sgc in 
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terms of independent fields. The Gaussian terms will define the bare diffuson and Cooperon 
propagators, while higher order terms will provide the corrections due to quantum inter- 
ference. To this end it is useful writing the unitary transformation U, introduced in Eq. 
Q, in an exponential form U = e^^^. The symmetry properties of W are the following: 
(i) anti-hermiticity W = —W\ (ii) invariance for charge conjugation CW'^C'^ = W\ (iii) 
non commutation with Qsp (i.e. with S3). Because of this last property, W has the energy 
structure W oc asi + [3s2 and Q can be written as 

Q = Qsp e"^ (12) 

To the above properties of W, one has to add a further constraint when considering the 
superconducting phase. Indeed, the presence in Eq. ((Tj) of the Cooper pairing term, which 
is proportional to T2S1 and determines the non conservation of the electric charge, requires 
the invariance condition [U, T2S1] = 0. This leads to 

[W,T2S,]=0 (13) 

thus reducing the number of massless modes in the superconductor. Indeed, in this case, the 
W field will be forced to be proportional in the energy space either to si or to S2. This sym- 
metry reduction plays a crucial role in the evaluation of the dephasing time, since it implies 
the vanishing of the singlet term of the electron-electron interaction in the superconducting 
systent^ii^, as explicitly discussed below. 

The final outcome for W fields both in metallic and in superconducting phase is summa- 
rized in the table of the appendix taken from Ref. 0. The table states whether the various 
components of the W-fields are real or imaginary and symmetric or antisymmetric. 

By inserting Eq.lfT^ in Eq.ffTIH). in which QS3 — > Xm^m and the traces are extended to 
many Matsubara frequencies as required by analysis of interacting case, and by taking into 
account the symmetries of the W^s the two-particle propagators read^ 

(^|'.,n^(q)^!.,.,(q')) = (±)^i^^5(q + q')^n^(9) 

(14) 

X {5Z6t% [±] <e + (-1)* [ s:%6'^_^ [±] 

where 5", T refers to the singlet-triplet components (i.e. W = Ws<Jo + i Wt-(^), i = 0, 1, 2, 3 is 
the index of the r matrices representing the particle-hole component (i.e., Ws/t = Wg^j,To + 
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i Ws/t-'t), abed and nmrq are replica and energy indexes respectively, and A„ = sign(c<j„). On 
the right side (±) applies depending whether the W-components are real or imaginary, and 
[±] whether they are symmetric or antisymmetric. Finally Dnm is the (Gaussian) propagator 



47riVo Dq^ + \uJn - UJm\ 
-D{q,UJn - UJm) 



with the diffusion coefficient D = as/{2No). 

The coefficients in Dnm will be corrected both by higher order quantum interference effects 
and by the electron-electron interactions. Quite generally, both these effects determine a 
renormalization of the diffusion coefficient. In particular, interaction effects result into a 
self-energy S(g,i7) which gives rise both to a renormalization of the diffusion coefficient, 
D ^ D, and to a renormalization of the frequencies by a factor Z. More interestingly 
the interactions determine also the appearing in D{q,Q) of a mass term (a non vanishing 
S(g = 0,i7 = 0)), which represents the inverse dephasing time r~^. When all corrections 
are included the renormalized propagator will then read 

4vr7Vo Dg2 ^ z \uJn - uJm\ + r- 



where D and Z indicate the renormalizations evaluated in Refs. while ^ will be 

explicitly computed in the following section. 



III. EFFECTS OF THE INTERACTION: CALCULATION OF THE DEPHASING 
TIME 

In the metallic phase the relevant interactions in the presence of disorder are the zero 
harmonic singlet and triplet amplitudes in the particle-hole channel Tg and Tt and the singlet 
s-wave amplitude in the Cooper particle-particle channel. In the d-wave superconducting 
phase the relevant residual Cooper interaction comes from the is-wave channel. In the 
following we will not consider this Cooper term which has been already analyzed in Ref. 

Indeed in the case of an additional repulsive zs-wave Cooper interaction, it can be shown 
with straightforward but lengthy calculation that its contribution to l/r^p is subleading with 
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respect to the triplet particle-hole contribution (Eq. below). This is because a repulsive 
coupling in the Cooper channel scales to zero and no relevant corrections survive from the 
particle-particle interaction. The situation would be completely different in the case of an 
attraction in the is-wave Cooper channel, leading to a. d + is instability, but the analysis of 
this regime, considered in Ref. Iisj, is beyond the scope of this paper. 

Starting from and Ft, and using the same formal steps for the metallic phase^^, the 
following additional contribution to Ssc is obtained^ 

S^nt[Q] = ^ E E (17) 

n,m,uj,a 1=0,3 

X {F, Tr(Q- +,r,ao)Tr(Q-+,,^r,ao) 
+F, Tr (Q-^^r, a)Tr(Q-^^,„r, a)} 

Notice that the Q matrices are diagonal in the replica space, since the interactions are 
present at each fixed disorder configuration. 

At this point we can make explicit the vanishing of the singlet contribution in Eq. (jl7p 
due to the additional constraint (jl3j) related to the absence of particle conservation in the 
superconducting phase. Indeed, as discussed in the previous section, the presence of the 
superconducting gap causes a reduction of diffusive degrees of freedom with respect to the 
normal case. More precisely, the charge-conjugation invariance, which stands both in normal 
and in superconducting systems, implies 



QTi,nm ~ ^ QTi,mn (l^) 

expliciting the charge- conjugation invariance for the singlet (S) and the triplet (T) particle- 
hole channels {i = 0,3), while in the superconducting case the following constraint also 
holds 



T2S1QT2S1 — -Q =^ Q'sLnm — i~^y^^ Qfi,-n-m 



QTi,nm — i ^y^^ Qri-n-m (l^) 



which follows from Eq.(|T^. 
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In general, if we consider the following transformation for the Q's 



Qnm ^ Q —m—n (^^) 

which is achieved by both permuting and changing the sign of the energy indices, we can de- 
compose every Q field in terms of its symmetric {sym) and anti-symmetric (ant) components 
under the transformation (plj) . Such a decomposition is useful, since in the superconduct- 
ing case the singlet and the triplet components of the Q^^-fields have definite symmetry 
properties under the transformation ((201) • particular, from Eq. (fTHj) and (fT^ one finds 
that 

Q'^S,sym = Q'T,ant = (21) 

On the other hand, because of energy conservation, Eq. p7|) can be written only in terms 
of the symmetric components, by rearranging the energy indices in the energy sum. As a 
result, only the triplet survives in Eq. (fTTj) . 

We are now in position to evaluate the dephasing time, i.e. the mass of the propagators 
which is generated by the electron-electron interaction. In the metallic phase this mass term 
comes from a specific one loop self-energy diagram containing a finite contribution from a 
branch cut. We get the same result here. To this end, after expanding the Q-field appearing 
in Sint[Q] in term of the H^- fields, we consider the four-W-field term 

Sint{W) = -32^ X]n,,m, Sa,fe,c S«=0,3 Sk,k',q -'^t 

X Tr (A„,W„tn^.(k)W^^%,(-k + q)r,a) (22) 

X Tr (A„3M/-^,(k')l^r.n4(-k' - q)^/^) K-n,,n,-ns 

By performing the contractions of the 14^— fields as it is shown in Fig. ^ and using 
the Gaussian propagator defined in Eq. (fT^ . we obtain the following contribution to the 
self-energy 

Sn,m = 37rT(l — Xm^n) Tt (1 " ^nj^ Xm) D nj^rn (23) 

ni 
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where we relabeled mi = m, n2 = n. This contribution is the counterpart in the supercon- 



ducting phase of the contribution coming from diagram [h) in Ref. 
following we take Um > and a;„ < 



221 To be specific, in the 



"i n,+l 





m.k,+q\ ? lm,k+q 



FIG. 1: (a) Four field term coming from the perturbative expansion of Smt in W (Eq. ((22)); (b) 
Diagrammatic representation of the self-energy term of Eq. H23|) 

As in the metallic phase, the dynamic structure of the interaction must be included in 
Eq. lj^Hj) . Indeed the one-loop self-energy is first order in t (the disorder parameter) but its 
calculation can be carried out to infinite order in Ft- This is achieved simply by replacing 
the static amplitude Tt in Eq. p3|) with the dynamic amplitude Tt^uj), given by^^*^^^ 



Dq^ + 






Dq^ + (1 - 


2r,) 


Uj\ 



where uj = ujn — cUni- 

The self-energy S, therefore, becomes 

with VL = UJm — UJn> ^■ 

In writing the energy limits in Eq. (j25|) we have considered the partial cancellation dis- 



cussed in Ref. 
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leading to the upper cutoff uo < cj^ + fi, so that finally the frequency 
sum has to be performed in the interval uJn < oj < uJn + Apparently this contribution 
is proportional to f2 and should be taken as a contribution to Z. However in carrying out 
the sum by analytically continuing on the real frequencies uj —iuj a branch cut at = 
appears. Along the branch cut Tt{p.,oj) and D{p,uj + Q) have different analytic continua- 
tions. Because of this, the sum does not vanish in the limit = and, in addition to the 
contribution to Z, the following mass term is also obtained 
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(26) 



ttNq sinh Pu 

X j j^;^^Im{rt{uj,p)}RDR{uj,p). 

Notice that the integral in Eq. is infrared divergent and takes contribution from 

momentum-energy scale < Dq^^uj < T. As pointed out by Fukuyama^^ '^/'t^ should 
enter selfconsistlently the diffusion propagator in the right side of Eq. ()26|). In practice, 
following Fukuyama we evaluate the integral with a low frequency cutoff of the order of 
Eventually we get: 



27r tT 



;i-r,) * 



In-^ (27) 



where 5 = 3tr^/(l — Ft) represents the cut-off coming from the self-consistency within 
logarithmic accuracy. 

This result indicates the existence of a direct proportionality between the energy scale 
and the temperature, with a coefficient of proportionality that contains the square of the 
scattering amplitude Tf. Comparing Eq. (j^7j) with the results for normal metals^i, we find 
that, despite the more involved Nambu formalism and the reduction of symmetry introduced 
by superconductivity, the triplet contribution to 1 /t^ in the d— wave superconducting phase 
has the same expression as in the normal phase. Conversely, as discussed above, one finds 
the vanishing of the singlet contribution, which is usually the most relevant contribution 
in the normal phase, due to the symmetries of the W fields in the superconducting phase. 
Notice also that usually other inelastic processes (i.e. electron-phonon or electron-electron 
interactions within the standard "clean" Fermi- liquid picture) lead to oc with p > 1 
so that they should be subleading at low temperature with respect to P7j) . 



IV. SIZE OF THE LOCALIZATION CORRECTIONS AND COMPARISON WITH 
EXPERIMENTS IN CUPRATES 

In this section we use the expression ()27|) to estimate the ratio t^/tq in the supercon- 
ducting phase of the cuprates at low temperature. To this end we express the parameters 
appearing in Eq. ()27|) in term of experimentally accessible quantities, within few reasonable 
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assumptions. The parameter Tt can be expressed as 



Tt = N^Vt (28) 

where Vt is the dimensional static scattering amphtude and for clarity we added the index 
s to the density of states of the superconductor, A^^q. An estimation of Vt is then obtained 
considering that its value should be essentially unaffected by the superconducting transition 
and that, assuming a weak- or intermediate-coupling for the metallic phase, one has Vt ~ 
1, i.e., Vt ~ where A*"^ is the metallic density of states. As a consequence, being 

A'q ~ l/irc'^vpVATo) « A"™ ~ kp/{2TiVF), the parameter Tt will read 

|r,| ^ i!^L ^ ^ ^ « 1 (29) 

where A ~ kpVA, being kp the Fermi momentum. 

Inserting Eq. pHj) and t = VpVAf^ivp + v^) ~ v^/i'^vp) in Eq. (j2I|) one gets 

TckpVp '^^F/A \ ron\ 

— - i:^^{kpVATo) - - — (Aro). (30) 

where we have dropped all the logarithmic corrections and Sp ^ kpVp. Notice that Eq. 

should be considered as an upper bound for Vt and consequently Eq. (jHn|l represents a 
lower bound for t^/tq. By comparing with the normal phase, where the largest contribution 
to is given by the singlet interaction^- and ^ utT ~ T / {kpVpTo), we find that 
t^/tq is enhanced in the superconducting phase by a factor Atq, which can be quite large if 
disorder is weak. 

At optimal doping the values for kp, vp and can be extracted by the ARPES data, 
while the elastic scattering time tq can be deduced from optical conductivity measure- 
ments. One gets kp ~ 0.8A~^, vp ~ 1.6eVA, vp/v^ ^ 19 (BSCCO) and vp/vA ^ 14 
(YBCO(123))^^»i^^. The evaluation of tq is the most problematic one. A rough estimate 
of 1/tq is obtained by extrapolating to lower temperatures l/r(T) derived from microwave 
measurements of the charge conductivity cr{T) (in the temperature range of 5 30 K) and 
taking I/tq = hniT^o^/riT). This leads to I/tq ~ OAK (YBCO(123)) and I/tq ~ 8K 
(BSCCO)M^. 

From these values we obtain the following estimates for the dephasing time in optimally 
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doped YBCO(123) and BSCCO 



4- lO^- 1 



YBCO(123) 




BSCCO 



(31) 



We can turn now to estimate the localization corrections to the transport coefficients 
according to Eq. 0: 



As stated before, the log factor contains the elastic scattering time. As a consequence, 
the localization corrections tend -rather surprisingly- to increase with decreasing disorder. 
The outcome of our calculation is that for both materials we would predict that localization 
effects should be visible in the temperature range we are interested in, i.e., hundreds of 
mK. In particular for a temperature of 1 K we get for 6as/crs = ^i^res/ i^res the values —0.32 
and -0.18 in the case of YBCO(123) and BSCCO. At 100 mK we get -0.36 and -0.23 
respectively. More than these specific values we want to stress that our calculation (which 
possibly underestimates the effect) gives a relative variation of as and Kres of about (or 
even more than) 20 ^ 30 per cent for temperatures lower than IK. The point is that the 
t = 'V/\/{2vf) is small (but not dramatically small) while t^^/tq is large enough to provide 
a sizeable correction. For YBCO(124) we expect similar values than YBCO(123) or even 
larger since these compounds are usually quite clean. 

These results are quite puzzling because, as we mentioned in the Introduction, no relevant 
localization effects are found in this temperature regime but for the case of YBCO(124) and 



The simplest possible explanation of such a discrepancy could be related to the existence 
in the real materials of other possible sources of inelastic scattering, beyond the diffusion 
enhanced electron-electron interaction (in particular scattering with some kind of quantum 
critical fluctuations^), and in principle this would be responsible of a further contribution 
to the dephasing effects. 

We can try to ge t the effective dephasing time from the temperature dependence of t(T) 
obtained in Ref. QandQ, as we did for tq. By assuming that the temperature dependence 
comes from inelastic processes and writing l/r(T) = I/tq + l/T:j„(T), we take as a rough 




(32) 



PCCO. 
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estimate for the dephasing time r^"^ ~ Tin(T), even though the two times are conceptually 
different. This procedure gives 



The comparison of Eqs. (jHHj) and ()3ip does not reveal any inconsistency in taking the 
diffusion enhanced electron-electron interaction as the predominant source of dephasing at 
low temperature in YBCO(123). In particular in the temperature range around IK our 
previous estimate of and the value of r^^^^ are of the same order of magnitude but display 
a different temperature dependence with the "experimental" estimate growing as in 
the low temperature regime. As a consequence the scattering due to the diffusion enhanced 
el-el interaction is dominant below T ~ 0.5 K and our previous theoretical estimate of Eq. 
()3H) has not to be corrected. On the contrary, in BSCCO the "experimental" displays 
the same 1/T behavior of the theoretical prediction, but with a much smaller coefficient, 
suggesting the presence of some other mechanism of inelastic scattering more effective than 
the diffusion enhanced electron-electron interaction. In this case it becomes more appropriate 
to use the "experimental" estimation of to evaluate the value of Sures/ i^res, the result being 
SKres/i^res — "0.03 (at IK) and -0.05 (at lOOmK). These values can be taken as an indication 
of negligible weak localization corrections in BSCCO because of strong decoherence effects. 

According to the above discussion and contrary to the experimental finding, significative 
localization effects would be expected in YBCO(123) as well as in YBCO(124), while the 
presence of a different source of strong inelastic scattering in BSCCO would explain the 
irrelevance of the localization corrections even at the lowest temperatures experimentally 
accessed. 

V. CONCLUSIONS 

In this paper we have dealt with the problem of the localization in high- Tc supercon- 
ductors. Their quasi 2-dimensional lattice structure would suggest that localization effects 
should become relevant at least in the low-temperature regime. However, as discussed in the 




BSCCO 



YBCO(123) 



(33) 
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Introduction, the low-temperature transport measurements in the cuprates depend strongly 
on the different materials and point out to the inadequacy of the prediction of a complete 
localization with sizeable weak localization precursor effects. Of course the electronic in- 
teractions could play in principle a relevant role in the explanation of the real material 
properties. 

A main effect of the interactions, which are a source of inelastic scattering, is to introduce 
a low-energy scale, the inverse of the dephasing time, which limits the quantum interference 
processes. Interactions also generate completely new corrections to transport, which could 
in principle compete with the localization effects. However they are expected to be negligible 
in the temperature-range of interest if disorder is weak^. 

We have derived the expression for the dephasing time in rf-wave superconductors 
within a non-linear a-model formulation of the disorder problem in the presence of inter- 
actions. Our result differs from the result for the metallic phase because of the vanishing 
of the singlet contribution, due to the symmetries of the superconducting phase. The main 
contribution to comes from the interaction in the triplet channel, while a repulsive 
residual interaction in the Cooper channel only produces a subleading term. 

The scale allows to estimate the size of the localization corrections to the heat (and 
spin) transport. Indeed we have estimated the localization corrections that are theoretically 
expected in the various materials. We have also considered the possibility that other sources 
of inelastic scattering could provide smaller dephasing times in real materials, and attempted 
to provide rough experimental estimates of l/r<^ to make a direct comparison with the data 
for thermal transport. The outcome of this comparison is that the estimated dephasing 
time is not small enough to reconcile the theoretical and the experimental findings in the 
cuprates. This is particularly evident from the failure of the theoretical predictions for 
the localization corrections in the case of YBCO(123). Moreover, as pointed out in the 
Introduction, the scenario of the low-temperature transport in cuprates was made even 
more involved by very recent experiments in LSCO^ and in YBCO(123)^, which claim for 
a strong doping- dependence of the thermal conductivity. 

The permanence of inconsistencies with the prediction of the localization theory, even 
when the electron-electron interactions and dephasing are taken into account, provides a 
clue of the relevance of some other effects which have been neglected in our analysis. For 
instance our theoretical predictions could be changed substantially by the proximity to 
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a d + is instability or to a Quantum Critical Point (which would enhance considerably 
the estimate of the energy scale by proximity to a nesting condition, by small q 

scattering and domain wall scattering, by the removal of the hypothesis of weak-coupling 
interactions, which could become questionable especially for the underdoped cuprates and 
large disorder. All the above effects will lower the temperature scale below which localization 
become sizeable. 

Finally we would like briefly to mention the relevant problem of magnetic impurities in the 
localization process. Since the scattering with magnetic impurities determines^*^ a vanishing 
of the localization corrections, a measurement of the low temperature transport coefficients 
in presence of magnetic impurities would be of a great relevance for understanding the role 
of the Anderson localization in the low-temperature transport of the cuprates, especially 
in the case of YBCO(124), where localization could be at work in providing a vanishing 
thermal transport. In YBCO(124) it would be also useful a systematic and careful analysis 
of magnetoconductance and thermal transport at small fields to compare with the findings 
in the cuprates which do not show localization. For these latter systems it would instead 
be important to assess the absence of magnetic impurities and local moments by magnetic 
measurements. 
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APPENDIX 

We report here explicitly the symmetries of the transverse massless modes W: 
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where n and m are odd integers which label the Matsubara frequencies c<j„ = Tirn and 
Um = Tnm, which have opposite sign {uJnOJm < 0). Notice that because of Eq. (fT!^ in the 
superconducting phase Wnm and W-n-m are not independent. 
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